Approach to the semiconductor cavity QED in high-Q regimes with q-deformed boson 
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The high density Frenkel exciton which interacts with a single mode microcavity field is dealed 
with in the framework of the q-deformed boson. It is shown that the q-defomation of bosonic 
commutation relations is satisfied naturally by the exciton operators when the low density limit is 
deviated. An analytical expression of the physical spectrum for the exciton is given by using of the 
dressed states of the cavity field and the exciton. We also give the numerical study and compare 
the theoretical results with the experimental results 



PACS number(s): 42.50 Fx, 71.35-y 



I. INTRODUCTION 

It's well known that the exciton system is a quasi- 
particle system. At the low density, excitons are approxi- 
mately treated as bosons which obey Bose statistics [Q-D . 
But when the density of the excitons become higher, the 
excitons, which somewhat deviate ideal bosons, are no 
longer ideal bosons. There are two ways of dealing with 
this problem : one way is to put these deviations into 
the effective interaction between the hypothetical ideal 
bosons and the exciton operators are still presented by 
the bosonic operators |4|j|]. Another way is the imple- 
ment of the atomic operators model 1^-^ . The question 
naturelly arises whether the exciton system is equivalent 
to the atomic system. 

In this paper we show that the exciton system could 
be described by q— deformed boson system which inter- 
polates between Bose system and Fermi system Q ; and 
the deformation parameter q is well defined by the total 
atomic particle number N^vaXhev than it is phenomeno- 
logical in the previous discussions .The concept of the q- 
deformed boson was even extensively applied due to the 
q-deformed boson realization of quantum group theory 
by different authors ten years ago pC|-p^. Since then, 
many physicists make efforts to find its real physical re- 
alizations. For example, they give some phenomenologi- 
cal investigations to fit the deformed spectra of rotation 
and oscillation for molecules and nuclei |p^-|l6|]. In our 
opinion, those investigations can be regarded as merely 
phenomenological because a g-deformed structure is pos- 
tulated in advance without giving it a microscopic mech- 
anism. In this paper it will be shown that a physical and 
natural realization of the (/-deformed boson is provided 
by the exciton operators, which was presented recently by 
Gardiner for the description of Bose-Einstein conden- 
sation (BEG). In fact, the similar quasiparticles scheme 
for particle-number conservation has already been intro- 



duced by Girardeau and Arnowitt almost 40 years ago 
|]l8| . The relationship between Gardiner's phonons and 
these quasiparticles has been discussed in a recent com- 
ment . 

Following these ways, we find that the exciton system 
also could be described by the q-deformed boson. When 
the density of the exciton (the particles excited in upper 
state) is low enough, it return to the ideal boson case. 
Using this theory, we could give a good explanation on 
the semiconductor cavity QED in high-Q regimes. What 
will be investigated here is the case that the total atomic 
particle number N is very large but not infinite. That 
is, we shall consider the effects of order o{l/N). And we 
shall focus on an algebraic method, a q-deformed boson 
algebra, of treating the effects of finite particle number. 
As it turns out, the commutation relations for the exciton 
operators will no longer obey the commutation relation of 
the Heisenberg-Weyl algebra but the g-deformed bosonic 
commutation relation 
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where the deformation constant q depends on the total 
atomic particle number. 

This paper is organized as follows. In section 2 we 
firstly deduce the q-deformed commutation relation for 
the exciton in the high-Q cavity in case of the large but 
finite lattice molecule number N . In section 3, only 
keeping the first order term of ;^,we model the Frenkel 
excitons in a micro-cavity as the dressed q-deformed bo- 
son system. In section 4, the quantum approach for an- 
gular momentum is used to obtain the eigen- values and 
eigen-function of the system under the one order approx- 
imation. The stationary physical spectrum of the system 
is calculated in the section 5. Finally we summarize our 
results with some comments. 



II. Q-DEFORMED BOSONIC ALGEBRA FOR 
EXCITON 

Gardiner's starting point @to introduce the exciton 
operators is to consider a system of the weakly interact- 
ing Bose gas. Without losing generality, we consider a 
thin molecular crystal film containing N identical two- 
level molecules interacting resonantly with a single mode 
quantum field. The intermolecular interaction is ne- 
glected. We assume that all molecules have equivalent 
mode positions, so they have the same coupling constant 
K. By using Dick model ||2^, we could write the Hamilto- 
nian under the rotating wave approximation as following: 



H = hn{S^ + a+a) + hK{aS+ + a+S-), 



(2) 



where, a is annihilation operator of the quantum cavity 



field and 
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S+ 
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where, s^(n) = i(|e„ >< e„| - |g„ >< 5„|), s+(n) = 
Icn >< 5n| and S-{n) = \gn >< e„| are quasi spin op- 
erators of the nth molecule. |e„ > and |(/„ > are the 
excited state and the ground state of n'th molecule. 

Consider the second quantization of the above model. 
Let bl and denote the creation and annihilation oper- 
ators for the atoms in the excited state and &J and bg for 
the creation and annihilation operators of the atoms in 
the ground state. The simplified Hamiltonian in second 
quantization reads 

H = hQ{blbe - blbg + a+a) + hK[ablbg + H.c.]. (4) 

Note that the total atomic particle number N = bjfee + 
b^gbg is conserved. For convenience we define rj — 1/N for 
large particle number. 

In the thermodynamical limit N ^ oo, the Bogoliubov 
approximation pl|j2^ is usually applied, in which the lad- 
der operators &|7^g of the ground state are replaced by 
a c-number ^/Nc, where Nc is the average number of 
the ininital condensated atoms. As a result Hamiltonian 
Eq.(^ becomes a two-coupling harmonic oscillator sys- 
tem 



Hh = nn 



+ a+a) + hK^/Nllabl + H.c.]. 



(5) 



However, this apporoximation destroyes a symmetry of 
the Hamiltonian Eq.(^, i.e., the conservation of the total 
particle number is violated because of [TV, TJf,] ^ 0. 

To avoide this problem, the exciton operators are de- 
fined as: 



(6) 



according to Gardiner ||17f|. These operators act invari- 
antly on the subspace spanned by bases \N;n) = 
\N — n,n) (n = 0, 1, . . . , N), where Fock sates (m, n — 
0,1,2,...) 

|m,n) = ^i=e6t«|o) 
Vm!n! 

spann the Fock space H2b of a two mode boson. 

A straightforward calculation leads to the following 
commutation relation between the exciton operaotr and 
its Hermitian conjugate: 



Kb\]=^-j^h% = f{blb,-r,), 



(7) 



with f{x;f]) = v'l + 2(l-2a;)?7-h if — rj. Keeping only 
the lowest order of ry for a very large total particle num- 
ber, the commutator above becomes 



or 



[bg,bl]=l-2T^b% 



(8a) 



(8b) 



with q = 1 — 277. This is exactly a typical q-deformed 
commutation relation. As — > cxd or g — > 1, the usual 
commutation relation of Heisenberg-Weyl algebra is re- 
gained. 

In the above discussion about the phonon excitation, 
we have linearized commutator h = f{b^b; 77) so that a 
q-deformed commutation rule was obtained. Essentially 
this linearization establishes a physical realization of the 
q-deformed algebra. However, if the total particle num- 
ber A^ is not large enough, then h can not be approxi- 
mated by a linear function. From the commutation rela- 
tions between h and bg, 6 J 



[h, bg 



N 



(9) 



we see that the algebra of exciton operators is a rescaling 
of algebra su{2) with factor A^. 

III. THEORETICAL MODEL 

Based on the above analysis about the algebraic struc- 
ture of exciton operator, we consider the case of the low 
density of atoms in excited sate for the Hamiltonian (2). 

Since the second quantization forms of S+ and S- are 

5+ = blbg, S- = b^gbc, it is straighforward to prove that 
s s 

the collective operator and are approximately 
considered as the simple bosonic operators as A^ ^ 00. 
These collective operators are called exciton operators. 
But in case of the high density of molecules in excited 
state with finite A^, many molecules are in the excited 
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states, the bosonic approximation could no longer work 
well. The Hamiltonian (2) is rewritten as the effective 
Hamiltonian in the form of q-deformed boson: 

H ^ hn{a+a + b+bg) + hg{a+bq + b+a) (10) 
with g = \/Nk, bq and 5+ satisfy q-deformed relation: 

[bq,b+]^bqb+ ~qb+bq^l, (11) 

where, 



(12) 



Here, the deformation parameter q is determined by the 
lattice molecule number. So q is no longer phenomeno- 
logical. 

Up to the first order approximation, b'^(bq) could be 
expressed as following 



,+ .+ b+b+b 

bt = b+ ^ , 

« 2iV ' 



b.^b 



b+bb 
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(13) 
(14) 



in terms of the general bosonic operators b'^{b) .So the 
Hamitonian H is rewritten in form of perturbation 



H = Ho + H' 



where 



Hq = hVt{a+a + b+b) + hg{a+b + 6+a), 



(15) 



(16) 



H' = —{2nb+b+bb + gb+b+ab + a+b+bb). (17) 

It is clearly that H' is equivalent to the attractive exciton- 
exciton collisions due to the bi-exciton effect and de- 
creased exciton-photon coupling constants due to the 
phase phase filling effect |24| . 



IV. APPROXIMATE ANALYTICAL SOLUTIONS 

To solve the Schroedinger equation governed by 
eq.(15),we implement the quantum angular momentum 
theory H 

If we define the angular momentum operators 



then 



Jz - i(a+a-6+6), 
J+ = a^fo, J- = ab^ . 



Jx = -{a+b + ab+), 
Jy^^ia+b-ab+). 



(18) 



(19) 



We rewrite the Hamitonian (5) 

Ha = nriN + 2hgJ^ = hQN + 2hge-'^-^y J.e'^-'y . (20) 

In terms of a S0{3) rotation h{flN + 2gJx) by e'^'^y. 
Note that the excitation number operator N = a^a + b^b 
is a constant under the a 5*0(3) rotation and 



= J,2 + J,2 + J| = ^(^-1) 



N N 



(21) 



is the total angular momentum operator. In terms of the 
eigen states of and 

(n + V+"^(b+V"^ 

\jm>= )ll_J^y ==|0>, (22) 
where the eigen values of the and Jz are 



(23) 



the eigen functions and the eigen values of Hq are con- 
structed as 

= e"'^^" \im >, Ef^^ = nnN + 2hgm. (24) 

Up to the first order approximation, the eigen values 
of H are obtained as 

Ejm = E^ll+ < mj\\e'i''yH'e~'^-^y\jm >, (25) 

whith their corresponding eigen functions 



HL 



nk ,/,(0) 



^(0) ^(o)^}:;^- (26) 

where, n and k present the subscript {jm') and {jm). We 
calculate the matrix elements of the pretubation Hamil- 
tonian H': 

< m' jWe'i-^y H' e~''^'^y\jm > 



= — f^V (j + + - 1) 

X \/{3 - m + l)(j - TO + 2)(5m_2,m' 

+ ^r!v/(.?+TO + l)(j+m + 2) 
X \/ [j ~ m){j - m - l)5„i+2,m' 



+ ^(2f^ - .9)V(j-^)(j+™ + l) 

h 



+ — (2r! + .g)V(j-m)(j+TO + l) 

X (j + m)Sm+l,m' 



+ ^(2r! + .g)v/(j+m)(j-TO + l) 
X {j + m- l)(5„i_i,™' 
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+ -(2r! - .g)v/(i + "i)(j-^ + i) 

X (j - ■m)S„i-i^„i' 

+ J^i^ + 9){j + m){j + m- l)S,n,,n' 

+ ^ 9)ij - - TO - l)S,n,,n' 

+ - m^)6ra,m' (27) 

so the eigen values of H are 

Ej^ = hflAf + 2m% + j^^if - rn'^) 
+ -^{^ + 9)ij + m){j + m - I) 
+ ^{n-g){j-m){j~m-l) (28) 

In general, we could obtain all eigen functions of H under 
one order approximation by using eq.(24), eq.(26), and 
eq.(27). So the time evolution operator of the system is 
written as: 

^ ^ . 

U{t) = exp{-i—t) = X! X! '2^P(~*~^0IV'im >< tpjr 
j=Q m=-j 

(29) 



b+bb^ 



where cuir, 



-j-^ and ujnrn 



It's evident 



that j is determined only by the initial state |i >. So we 
have 



S{u;) = 



27 



72 + (cj - LOhn)' 



(33) 



Noting that we have passed the transient terms and 
slowly variation terms. This equation gives the station- 
ary physical spectrum in terms of the system eigenvalues 
and eigenstaes. If < m! j'\b^\jm >^ 0, then we have 
j' — j + and m' = TO — i. So the equation (25) is 
rewritten as: 



27 



72 + (tJ - UJhn}- 



b+b+b 
2N 



(34) 



V. FLUORESCENCE SPECTRUM OF THE 
EXCITON 

We firstly give an analytic expression for the physical 
spectrum of the q-deformed exciton in terms of the Fock 
state of the quantum field and the exciton. The standard 
definition of the physical spectrum is ||2^ 

5(a)) = 27/" dti f di2e-(^-^"'^(*-*'^e-(''+'")(*-*i)G(ti,i2; 
JO Jo 

(30) 

where 7 is the half-bandwith of spectrometer which is be- 
ing used to measure the spectrum, and t is time length of 
the excitation in the cavity. G{ti, 12) is dipole correlation 
function, and 

G(ti,t2) -< i\U+{t2)b+U{t2)U+{h)bqU{h)\i >, (31) 

where, \i > is any initial state of the system. By using 
the bosonic approximation and substituting eq.(29) into 
eq.(31), the dipole correlation function is rewritten as 
following: 

G{h,t2) l\U+{t2)b+U{t2)U+{h)bgU{h)\i > 

h+b+h h+bh 

= < t\U+{t2){b+ + ^^)U{t2)U+{h){b + -^))U{h)\z > 



The eigenvalues determine the position of the spec- 
tral component and | < > P| < ipji\{b^ + 
^-2^)IV'(j-i)m > P determine the intensity of the spec- 
tral lines. 

In terms of the experimental condition of the refer- 
ence the bare excitons could be prepared by reso- 
nant femtosecond pulse pumping. If we prepare the ini- 
tial state being in J\f = 1, then the eq.(34) shows that 
the emission spectrum of the JV = 1 to the Af = transi- 
tion has double peaks structure which is exactly equal to 
that of the two-level atomic system. When the pumping 
power is increased , the emission spectrum is quiet dif- 
ferent from the case of the two-level atomic system. For 
example, if the system initially is in Fock state Af = 2 , 
then we have: 



27 



<^ii\ib+ 



b+b+b 



2N 



I < > P 
)|^i,„ > p (35) 



From this equation we know when J\f = 2 there are six 
peaks in the emission spectrum which are expected. Al- 
though there are three different initial state, they have 
similar spectrum shape. As is shown in Fig.l, this sextet 
structure is contrast to the triplet structure in the emis- 
sion spectrum from strong pumped two- level system [p7[ . 
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VI. CONCLUSION 
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FIG. 1. f2=1562 mev, N=100, g=20 mev, 7=0.1 mev 

When the molecular number of the system is increased, 
Such as there are 10000 molecular in the system, the 
other conditions are the same as that of the Fig.l, then 
the coupling between the molecular and the cavity field 
is weak. There are two peaks in the emission spectrum 
(Fig. 2). In this case Bose approximation is good. 

q 




1500 1525 1550 1575 1600 1625 1650 



FIG. 2. fi=1562 mev, N=10000, g=20 mev, 7=0.1 mev 



It has been shown that the higher density Prenkel 
cxciton coupling to a single mode high-Q microcavity 
field can be described by the quantum dynamics for the 
dressed q-deformed boson. Keeping the first order term 
of jj ,thc high density Fronkcl cxciton naturally obycs the 
q-deformed commutation relation. Based on this obser- 
vation the quantum theory of angular momentum is used 
to obtain the cigcn- values and cigcn-fuuction of the cxci- 
ton system under the one order approximation. Compar- 
ing with the usual approach for Prenkel exciton dynam- 
ics our Hamiltonian is Hcrmitian and closed in form. An 
analytical expression for the stationaryphysical spectrum 
for the exciton is obtained by using of the dressed states 
of the cavity field and the exciton. 
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